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Abstract — In this paper, an evolutionary optimization 
technique, Particle Swarm Optimization with Constriction 
Factor and Inertia Weight Approach (PSOCFIWA) is adopted 
for the complex synthesis of three-ring Concentric Circular 
Antenna Arrays (CCAA) with non-isotropic elements and 
without and with central element feeding. It is shown that by 
selection of a fitness function which controls more than one 
parameter of the array pattern, and also by proper setting of 
weight factors in fitness function, one can achieve very good 
results. For each optimal design, optimal current excitation 
weights and optimal radii are determined having the objective 
of maximum Sidelobe Level (SLL) reduction. The extensive 
computational results show that the CCAA designs having 
central element feeding with non-isotropic elements yield 
much more reduction in SLL as compared to the same not 
having central element feeding. Moreover, the particular 
CCAA containing 4, 6 and 8 number of elements in three 
successive rings along with central element feeding yields 
grand minimum SLL (-46.4 dB). Standard Particle Swarm 
Optimization (PSO) is adopted to compare the results of the 
PSOCFIWA algorithm. 

Keywords — Concentric Circular Antenna Array; Non-isotropic 
Elements; Particle swarm optimization; Non-uniform 
Excitation; Sidelobe Level; First Null Beamwidth 

I. Introduction 

A Concentric Circular Antenna Array (CCAA) is an array 
which contains many concentric circular rings of different 
radii and a number of antenna elements located on the 
circumference of each ring [1, 2]. CCAA has received 
considerable interest due to its symmetricity and 
compactness in structure. Since a concentric circular array 
does not have edge elements, directional patterns synthesized 
with a concentric circular array can be electronically rotated 
in the plane of the array without a significant change of the 
beam shape. CCAA offers great flexibility in array pattern 
synthesis and design both in narrowband and broadband 
applications. These salient features of CCAA have made it 
indispensable in mobile and communication applications. 
These antenna arrays have improved the performance of 
mobile and wireless communication systems through efficient 
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spectrum utilization, enhancing channel capacity, extending 
coverage area and tailoring beam shape etc. However, antenna 
array design can not be done on an arbitrary basis as this 
may lead to pollution of the electromagnetic environment 
and more importantly, wastage of precious power. The later 
may prove fatal for power-limited battery-driven wireless 
devices. This has encouraged a lot of research work [3-11] in 
the field of optimisation of antenna structures, all having a 
common objective of bridging the gap between desired 
radiation patterns with what is practically achievable. 

The goal of this paper is to optimize current excitation 
weights and radii in order to minimize the SLL, hence working 
towards the improvement of antenna array pattern. In this 
approach the structural geometry of the antenna array is 
manipulated while preserving the gain of the main beam. This 
paper presents the design of concentric circular antenna array 
with non-isotropic elements. CCAA designs having non- 
isotropic elements (i) with and (ii) without central element 
feeding [12-13] have been considered in this paper. Many 
synthesis methods are concerned with suppressing the SLL, 
as the shape of the desired pattern can vary widely depending 
on the application. In this paper an optimal CCAA is designed 
with the help of Standard Particle Swarm Optimization (PSO) 
[14] and Particle Swarm Optimization with Constriction Factor 
and Inertia Weight Approach (PSOCFIWA) [15-16] 
techniques. The array factors due to optimal radii and non- 
uniform excitations of non-isotropic elements in various 
CCAA designs are examined to find the best possible design 
set. 

II. Design Equation 

Fig. 1 is an illustration of the general configuration of 
CCAA having M concentric circular rings, where the m th (m = 
1,2, . . ., M) ring has a radius r m and the corresponding number 
of elements is N m . If all the elements (in all the rings) are 
assumed to be isotropic sources, the radiation pattern of this 
array can be written in terms of its array factor only. But the 
elements considered in this paper are non-isotropic. When 
the actual elements are non-isotropic, the total field can be 
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formed by multiplying the array factor of the isotropic sources 
with the field of a single element. 

The far-field radiation pattern of the array ( FF(<j>, I, r) ) is 
then given by [1] 

FF (0,7, r) = EP (<f>,I,r).AF (<j>,I,r) (1) 
where EP (</>,! ,r) gives the element pattern (in this work; 
the elemental pattern is considered to be cos(0) type) while 
AF(<f), I, r) is the array factor of the isotropic source. 

With reference to Fig. 1, the array factor, AF{fff,I,r) for the 
CCA A in x-y plane maybe written as [12-13]: 

AF (<j>,I ,r) = X Z ex P[X^„ sin # cos ~<l>mi) + «™ ■)] (2) 
where 

I mj = current excitation of the i th element of the m th ring; 
k = In I X ; X being the signal wave-length; 

<f) mi =angular separation between elements measured from 
positive x-axis. 

If the elevation angle, Q = 90° then (2) may be written as a 
periodic function of (j) with a period of 2n radian. The angle 
is the element to element angular separation measured from 
the positive x-axis. The elements in each ring are assumed to 
be uniformly distributed, and are also obtained from [12] as: 

(j> mi =2n{{i-\)l N m ) (3) 

a mi =-Kr m cos^) -^ mi ) (4) 

<f> is the value of (j) where peak of the main lobe is obtained. 
After defining the far-field radiation pattern, the next step in 
the design process is to formulate the objective function 
which is to be minimized. 

The objective function "Cost-Function" (CF) may be 
written as : 

\FF(<j> msll ,I mi ,r m ) + FF{<j> mi ., 2 ,I mi ,r m 



CF =W F1 x 



+ 



W F2 x (fnbw 



\FF(^,I ml ,rJ 

FNBW (/„, = 1)) 



(5) 



computed 

FNBW is the First Null Beam Width. W Fl andW F2 are the 

weighing factors. Minimization of CF means maximum 
reductions of SLL both in lower and upper bands and lesser 
FNBW computtd as compared to FNBW (/„„ = l). The 
evolutionary techniques described below are employed for 
the optimisation of the current excitation weights as well as 
the radii of the concentric rings. This results in the 
minimization of and hence reductions in both SLL and FNBW. 
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Figure 1. Concentric circular antenna array (CCAA) 

III. Evolutionary Techniques Employed 

A. Particle Swarm Optimization (PSO) 

PSO is a flexible, robust population-based stochastic 
search/optimization technique with implicit parallelism, which 
can easily handle with non-differential objective functions, 
unlike traditional optimization methods. PSO is less 
susceptible to getting trapped on local optima unlike GA, 
Simulated Annealing, etc. Eberhart and Shi [14] developed 
PSO concept similar to the behavior of a swarm of birds. PSO 
is developed through simulation of bird flocking in 
multidimensional space. Bird flocking optimizes a certain 
objective function. Each particle (bird) knows its best value 
so far (pbest). This information corresponds to personal 
experiences of each particle. Moreover, each particle knows 
the best value so far in the group (gbest) among pbests. 
Namely, each particle tries to modify its position using the 
following information: 

• The distance between the current position and the pbest 

• The distance between the current position and the gbest 
Similar to GA, in PSO techniques also, real-coded particle 

vectors of population n p are assumed. Each particle vector 
consists of components or sub-strings as required number 
of current excitation weights and three radii, depending on 
the number of current excitation elements in each CCAA 
design. 

Mathematically, velocities of the particle vectors are 
modified according to the following equation [12-13]: 

Vj- k+1) = w * Vj- k) + C, * rand, * {pbestf ] - S\ k) ) 

+ C 2 * rand 2 *(gbest (k) - S\ k) ) ' V " 

where vf- 1 ' is the velocity of i th particle at k th iteration; w is the 
weighting function; C, and C 2 are the positive weighting fac- 
tors; rand, and rand 2 are the random numbers between 
and 1; sj k ^ is the current position of i th particle vector at k th 
iteration; pbestf" 1 is the personal best of i th particle vector at 
k* iteration; gbest^ is the group best of the group at k th 
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iteration. The searching point in the solution space may be 
modified by the following equation: 

5 = $(*) + yM (7) 

The first term of (6) is the previous velocity of the particle 
vector. The second and third terms are used to change the 
velocity of the particle. Without the second and third terms, 
the particle will keep on "flying" in the same direction until it 
hits the boundary. Namely, it corresponds to a kind of inertia 
represented by the inertia constant, w and tries to explore 
new areas. 

B. Particle Swarm Optimization with Constriction Factor 
and Inertia Weight Approach (PSOCFIWA) 

Steps of PSOCFIWA as implemented for optimization of 
current excitation weights and location of elements are 
adopted from [15-16]. 

IV. Simulation Results and Discussions 

This section gives the computational results for various 
CCA A designs obtained by PSO and PSOCFIWA techniques. 
For each optimization technique ten three-ring (M=3) CCAA 
designs are assumed. Each CCAA maintains a fixed optimal 
inter-element spacing between the elements in each ring. The 
limits of the radius of a particular ring of CCAA are decided 
by the product of number of elements in the ring and the 
inequality constraint for the inter-element spacing, 

d,(d e[A/2,l]) . For all sets of experiments, the number of 

elements for the inner most ring is N , for outermost ring is 
N,, whereas the middle ring consists of N, number of elements. 

For all the cases, O = 0° is considered so that the centre of 
the main lobe in radiation patterns of CCAA starts from the 
origin. Since PSO techniques are sometimes quite sensitive 
to certain parameters, the parameters should be carefully 
chosen. Best chosen maximum population pool size, n p =120, 
maximum iteration cycles for optimization, N =400. The PSO 
and PSOCFIWA generate a set of optimized non-uniform 
current excitation weights and optimal radii for each design 

set of CCAA. I nn = 1 corresponds to uniform current excitation. 
Sets of three-ring CCAA (N , N,, N.) designs considered for 
both without and with central element feeding are (2,4,6), 
(3,5,7), (4,6,8), (5,7,9), (6,8,10), (7,9,11), (8,10,12), (9,11,13), 
(10,1 2, 13) and (11,13,15). Some of the optimal results are shown 
in Tables II- V. Table I depicts SLL values and FNBW values 
for all corresponding uniformly excited (=1) CCAA design 
sets. 

Analysis of Radiation Patterns of CCAA 

Figs. 2-3 depict the substantial reductions in SLL with 
optimal non-uniform current excitation weights and radii, as 
compared to the case of non-optimal uniform current 
excitation weights and radii (considering fixed inter-element 
spacing, d=X/2) using PSO and PSOCFIWA optimization 
techniques respectively. All CCAA design sets having central 
element feeding (Case (b)) yield much more reductions in 
SLL as compared to the same not having central element 
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feeding (Case (a)). It is depicted from Tables II-III, SLLreduces 
to -34.8 dB for case (a) and -42.8 dB for Case (b) using standard 
PSO algorithm for the CCAA having N =4, N,=6, N 3 =8 
elements (Set No. III). Further for Tables IV- V, SLLreduces to 
-40.63 dB for Case (a) and -46.4 dB (grand minimum SLL as 
determined by PSOCFIWA for Case (b) and non-isotropic 
elements. This optimal design set along with central element 
feeding yields grand maximum SLL reduction among all the 
sets. The above reductions of SLL are determined by 
comparing Table I with Tables II- V. From comparison between 
Table II with Table IV and Table III with Table V it is clear that 
PSOCFIWA consistently yields better optimal results than 
standard PSO. 

Table I. SLL And FNBW For Uniformly Excited ( I mj =1 ) CCAA 
Sets 



Set 
Nd. 


No. ot 
elements in 
each rinES 

(Nl.N^Nj,) 


Without central 
element (Case (a)) 


Wth central 
element (Case (h}) 


SLL 


FNBW 
(dee) 


SLL 
(dB) 


FNBW 
(dad 


I 


2,M 


-14. B2 


12B.52 


-2C.J5 


140.0 


n 


3 5 7 


- J .:\44 


107.32 




116.5 


in 


4. 6. S 


-13.97 


rj.'A 


-12 31 


95.3 


IV 


5. 7. 9 


-16.2b 1 


7B.1 


-15.59 


SI. 54 


v 


6. S. 10 


-14.44 


6B.36 


-16.12 


71.22 


VI 


7. 9, 11 


-:<.:: 


60.92 


-14.56 


63.2 


VII 


S. 10. 12 


-12 .15 


54.82 


-13.37 


56.34 




9. n. 13 


-1142 


50.02 


-12.49 


51. IB 


IX 


10, 12. 14 


-10. B6 


46.02 


-11. B2 


47.15 


X 


11. 13. 1? 


-10.42 


4:.:: 




43.16 



Table II. Current Excitation Weights, Radii, SLL and FNBW For Non- 
Uniformly Excited CCAA SETS (Case (a)) Using Standard PSO 



"ETT 
No. 



TIT 



C • %i a -. — » }j 
(ri, ii. rj) in 1 



0.7279 0.4510 
0.3637 0.3369 
3S12 0.1436 
0.2373 0.1150 
0.1B57 



3.6SSI 0.4139 
0.1761 0.321B 
0.2199 0.4047 
0.1933 0.3614 
0.1432; 



0.32: 

3926 
0.4137 
0.7292 
0.4*75 
1711 
0.0093 



B 5529 
0.E2B6 
0.6174 

;.:::4 

0.2144 
0.6595 
0.0000 



1 0164 
3.B61- 0.B25O 
0.4013 0.0717 
0.602S 
0.6373 0.0000 
0.O49S 0.3955 
0.4130 0.2356:. 



0.4521 0.B070 1.1075 
0.741B 0.015B 0.6156 
0.5E43 0.2871 0.5113 0,61 3E 
0.7486 0.0001 0.7319 
0.4079 0.4306 0.12 S3 0.1735 
0.2601 0.3B91 0.0153 0.4770 
0.7203 0.5101 0.0036 0.3936 
0.470 5 0.1752 0.7B71 0.0727 
0.4921 0.4536: 

0.7021 0.9603 1.4334 



~STJT 

(dB) 



T34T" 



FNBW 
(dei) 



"B4T06" 



75.7B 



53^F 
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Table III. Current Excitation Weights, Radii, SLL and FNBW For 
Non-Uniformly Excited CCAA SETS (Case (b)) Using Standard PSO 



"TiT 
No. 



TIT 



TIT 



0.972(5 
VJUU 
0.8036 
0.19:6 
0.2292 
0.3446 



(ft- rj. rj) in. % 



.0= 



0.B091 

0.7994 

0.3990 

0.3245 

0.5475 



O?07 
0.8063 
0.2073 


0:. 
0.8561 



.1988 
.394 B 
3241 



0.1019 
0.1133 
0.2 
0.1 

0.24B7 
0.2 
0.157B 
0.5051 



0.1320 


594 0.4172 
8*4 0.290* 
0.2929 
672 0.22 5 5 



1*32 
2*90 
0292 
01ZS 
0315 



0.4405 
0.6409 
0.5004 
0.2*86 
0.1744 
0.0401 
0.5299 
0.0009 
0.70*3 
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0.0179 



0.165B 
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0.625* 

0.4*20 

1.0696 



1 2015 
0.3797 0.1236 
0.5578 
0.1746 0.30 87 



0.214* 

0.4427 

0.0B33 

0.0445 
0.3B4B 
1.3785 



0.2409 
0.1972 
0.3*30 
0.4*37 
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Figure 2. Radiation patterns for a uniformly excited CCAA and 
corresponding standard PSO and PSOCFIWA based non-uniformly 
excited CCAA Set No. Ill, Case (a) 
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Figure 3. Radiation patterns for a uniformly excited CCAA and 
corresponding standard PSO and PSOCFIWA based non-uniformly 
excited CCAA Set No. Ill, Case (b) 



Table IV. Current Excitation Weights, Radii, SLL and FNBW For Non- 
Uniformly Excited CCAA Sets (Case (a)) Using PSOCFIWA 



"5iT 
No. 



Tir 



TIF 



(ri, rj, rj) in 1 



0.41B7 0.4324 0.4253 0.4641 
0.5040 0.4870 0.7803 0.4756 
0.4 SB* 0.7873 0.3825 0.5105 
0.3685 0.3432 0.4948 
0.363B 0: 
0.30B0 0.5393 0.9312 



0.7 824 
0.73*9 
0.4351 
0.5887 
0.4185 
0.4241 



0.B597 
0.0503 
0.2771 
0.1:56 
0.:19: 



0.744B 0.7E6: 
0.5666 
03951 
0.2367 0.4*97 
0.1200 0.1788 
0.274B 0.1896; 



0.5087 0.B147 1.2364 



0.3745 
0.3072 
0.3032 
0.58*6 
0.3824 
0.495* 
0.3770 



0.01:: 


0.1700 
0.4136 
0.7337 


0.1:51 
0.19:0 
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0.0821 
0.0393 
0.3504 
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0.4814 
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0.2573 
0.2053 
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0.4B0B 
0.2793 



0.7407 1.0453 1.3BB9 
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Table V. Current Excitation Weights, Radii, SLL and FNBW For Non- 
Uniformly Excited CCAA Sets (Case (b)) Using PSOCFIWA 



"5iT 
No. 



"nr 



0.3036 
0.34*1 
0.5187 
0.2292 
0.2381 



(ri, rj. rj) in a 



0.3017 0.3472 
0.52 58 0.5174 
0.5254 0.3787 
0.2593 
0.2520 0: 



0.3106 
0.3B49 
0.2472 
0.2573 



0.3313 0.5103 0.8743 



TJ 0.319B 0" 

0.0927 0.2215 0.4213 
0.6739 0.3907 0.2300 
0.6072 0.0790 0.2683 0.2745 
0.1580 0.1401 0.2375 
0.2903 0.0511 0.01B2; 
0.5311 0.7465 1.1001 



0.1302 
0.5841 
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0.5501 
0.0049 
0.2030 
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0.1214 
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Convergence profiles of PSO and PSOCFIWA 

The minimum CF values are plotted against the number 
of iteration cycles to get the convergence profiles for the 
optimization techniques. Figs. 4-5 show the convergence 
profiles for PSO and PSOCFIWA for Set No. Ill, Case (b) 
CCAA respectively. PSO yields suboptimal higher values of 

CF but PSOCFIWA yields true optimal (least) CF values 
consistently in all cases. With a view to the above fact, it may 
finally be inferred that the performance of PSOCFIWA 
technique is better as compared to standard PSO. All 
optimization programs are written in MATLAB 7.5 version 
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on core (TM) 2 duo processor, 3.00 GHz with 2 GB RAM. 

1 1 1 1 1 




50 100 150 200 250 300 

It ntio n Cyc le 

Figure 4. Convergence curve for standard PSO in case of non- 
uniformly excited CCAA (Set No. Ill, Case (b)) 

s | 1 1 1 1 1 1 




50 100 1S0 2 00 250 300 

Iteration Cycle 

Figure 5. Convergence curve for PSOCFIWA in case of non- 
uniformly excited CCAA (Set No. Ill, Case (b)) 



V. Conclusion 

This paper illustrates how to model the optimal design of 
non-uniformly excited CCAAs with optimal radii and non- 
isotropic elements for maximum SLL reduction using Standard 
PSO and PSOCFIWA. Computational results reveal that the 
optimal design of non-uniformly excited CCAA offers a 
considerable SLL reduction with respect to the corresponding 
uniformly excited and d=e/2 inter-element spacing CCAA. 
The main contribution of the paper is threefold: (i) All CCAA 
designs having central element feeding yield much more 
reductions in SLL as compared to the same not having central 
element feeding, (ii) The PSOCFIWA technique outperforms 
significantly as compared to standard PSO in terms of solution 
optimality. iii) The CCAA design having Nj=4, N 2 =6, N 3 =8 
elements along with central element feeding and non-isotropic 
elements gives the grand minimum SLL (-46.4 dB) as compared 
to all other design sets as determined by PSOCFIWA, which 
one is thus the grand optimal design set among all the three- 
ring designs. 
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